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The Bianchi IX mixmaster model is quantized in its non-diagonal form, imposing spatial diffeo- 
morphism, time reparametrization and Lorentz invariance as constraints on physical state vectors 
before gauge-fixing. The result turns out to be different from quantizing the diagonal model ob- 
tained by gauge-fixing already on the classical level. For the non-diagonal model a generalized 
9-dimensional Fourier transformation over a suitably chosen manifold connects the representations 
in metric variables and in Ashtekar variables. A space of five states in the metric representation is 
generated from the single physical Chern-Simons state in Ashtekar variables by choosing five dif- 
ferent integration manifolds, which cannot be deformed into each other. For the case of a positive 
cosmological constant A we extend our previous study of these five states for the diagonal Bianchi 
IX model to the non-diagonal case. It is shown that additional discrete (permutation) symme- 
tries of physical states arise in the quantization of the non-diagonal model, which are satisfied by 
two of the five states connected to the Chern-Simons state. These have the characteristics of a 
wormhole groundstate and a Hartle-Hawking 'no-boundary' state, respectively. We also exhibit a 
special gauge-fixing of the time reparametrization invariance of the quantized system and define an 
associated manifestly positive scalar product. Then the wormhole ground state is left as the only 
C^") , normalizable physical state connected to the Chern-Simons state. 
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I. INTRODUCTION 



> 

On 
O 

Quantum general relativity has advanced over the last decade at a remarkable and accelerating pace. The in- 
traduction of Ashtekar's new variables |j]-[3| to replace the earlier metric representation soon afterwards led to the 
discovery of a formally exact physical state of the canonically quantized theory with non- vanishing cosmological con- 
stant, the Chern-Simons state [§J^]. Then the introduction of the loop representation || permitted to reexpress the 
Chern-Simons state as a topological invariant of framed loops on 3-space [Q, the Kauffman bracket It also led 
to the discovery of further physical states lacking, however, one important general property of the Chern-Simons 
state, namely a well-defined non-degenerate space-time in the classical limit. The Chern-Simons state semiclassi- 
cally describes a deSitter (or anti-deSitter) space-time for positive (negative) cosmological constant, respectively [Q. 
Subsequent important advances were the introduction of spin-networks Q and quantum spin-networks [|l0| as a dis- 
cretizcd description of 3-space in which areas and volumes are quantized in Planck units and which furnish yet another 
representation of physical states. 

The choice of different variables and representations has therefore played, and continues to play, a crucial role in the 
development of the theory. It is not always clear, however, to which extent the different representations are equivalent 
to each other. This question is particularly relevant for the connection between the metric representation and the 



representation in Ashtekar variables. In fact it was shown in |11| that the two representations are, in general, not 
equivalent. One may therefore wonder: Has the Chern-Simons state, the only known physical state with a well-defined 
classical limit, a counterpart in the metric representation and is it unique? In general, the answer to this question is 
still unknown. 

Recently we examined this question for a spatially homogeneous mini-superspace model of Bianchi type IX with 
diagonal metric tensor in the case of positive |l^] and negative jl3| cosmological constant A. A diagonal form of 
the metric tensor can be assumed in the classical theory without restriction of generality, because it is a permissible 
gauge-fixing condition for the remnant of the spatial diffeomorphism group in the Bianchi type IX model. Taking 
also the matrix of Ashtekar variables as diagonal amounts in addition to a gauge-fixing of the Lorentz gauge group. 
Quantizing such a diagonal model therefore means to apply gauge-fixing of the diffeomorphism group and the Lorentz 
group before quantization. The result of our study of the diagonal Bianchi type IX model was that actually five distinct 
physical states in the metric representation are generated by transforming the Chern-Simons state from Ashtekar 
variables to metric variables. This change of representation takes the form of a generalized multidimensional Fourier 
transformation in the space spanned by the complex Ashtekar variables along arbitrary paths or, more precisely, 
integration manifolds with boundaries pinned by the condition that partial integration without boundary terms must 
be allowed. 
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In the present work we take our investigation of this basic question a step further and examine the metric repre- 
sentations of the Chern-Simons state for the non-diagonal Bianchi type IX model. Why is this step interesting, and 
why can't the results of our study of the diagonal case not just be taken over? Quantizing the non-diagonal model 
amounts to interchanging the steps of gauge-fixing and quantizing: the gauge-fixing of the spatial diffcomorphism 
group and Lorentz group is now done after the quantization. The result is, in general, not the same as before the 
exchange of these steps. In fact, in section III A we give a very simple example from quantum mechanics which 
shows that gauge-fixing after quantization is preferable, in general, because it takes all symmetries into account in 
the quantization process in a natural manner. The example also makes clear that the result, in general, differs from 
that of first gauge-fixing then quantizing by (i) quantum corrections in the Hamiltonian, (ii) weight-functions in the 
naturally defined scalar product and (iii) additional discrete symmetry requirements to be satisfied by the solutions. 

We shall see that all three differences also show up in our study of the non-diagonal Bianchi type IX model. In 
fact the noncommutativity of gauge-fixing and quantization in the Bianchi type IX model was previously discussed by 
Major and Smolin within the framework of path integral quantization |16 17]]. Our investigation here differs by the use 
of canonical quantization and, in particular, by studying a special quantum state, the Chern-Simons state. Comparing 
the two ways of quantization we find that on the leading scmiclassical level the solutions for the non-diagonal model 
are the same as in the diagonal case. However, in the next to leading semiclassical order they already differ due 
to the quantum corrections in the Hamiltonian and due to changes in the naturally defined scalar product. Even 
more importantly, the additional discrete symmetry requirements are only met by two of the five linearly independent 
states, leaving just two physical states in the metric representation of the non-diagonal model which are generated by 
the Chern-Simons state, a generalized wormhole state and a Hartle- Hawking state |T^,|l5|. It is remarkable that states 
of both kinds are related to the Chern-Simons state and that the two different semiclassical boundary conditions 
singling out either one of them can still be satisfied at this stage. 

In a final step we also gauge-fix (after quantization) the time reparametrization invariance and introduce a manifestly 
positive scalar product on the space of physical states. As will be shown in a separate paper only the generalized 
wormhole state is normalizable in this scalar product. 

The rest of this paper is organized as follows: In section II we define our notation and set up the metric representation 
of the constraint equations of the non-diagonal model (II A), extract a well-known exact solution for A = 0, the 
wormhole state (II B), and give the representation of the constraint equations in Ashtekar variables (II C). In section 
III we discuss in detail the differences between the diagonal and the non-diagonal Bianchi type IX model. We first give 
a simple quantum mechanical example, the harmonic oscillator in two dimensions with a rotational gauge symmetry 
(III A). Then the corresponding comparison for the Bianchi type IX model is given (III B). We also mention here 
briefly a gauge-fixing of the time reparametrization symmetry and an accompanying physical inner product on the 
space of physical states (III C) but a more detailed presentation of this point is beyond the scope of the present paper 
and will be given separately In section IV the transformation of the Chern-Simons state of the non-diagonal 

Bianchi type IX model from the Ashtekar representation (IV A) to the metric representation (IV B) is given, leading 
to a general integral representation over a suitably chosen 3-dimensional manifold S 3 . In section V various asymptotic 
forms of this integral representation are evaluated. The limits considered are first h — > (V A) then either A — > 
(V Al) or A a 2 — > (V A2), where a is the geometrical mean of the scale parameters, and A — > (V B) without 
necessarily taking a second limit. In section VI we exhibit for the general case, i.e. without taking asymptotic limits, 
five possible and topologically distinct choices of the integration manifolds S 3 leading to five exact solutions of the 
constraint equations in the metric representation. We also discuss their relation with the asymptotic results of section 
V and their normalizability with respect to the inner product of section III C. Our conclusions are then summarized in 
section VII. Three appendices deal with certain technical details. The limit A — > 0, the results of which are discussed 
in section V B, is quite subtle and therefore done in some detail (A). It leads to a nice integral representation of 
the vacuum (i.e. A = 0) solutions (B). Finally we check certain required nontrivial continuity and differentiability 
properties of the integrand of the integral representation of section IV B on the integration manifolds S 3 (C). 



II. THE QUANTUM CONSTRAINT EQUATIONS 

In this section we shall set up our notation and give a brief derivation of the quantum Einstein equations for the 
homogeneous Bianchi IX model. While the classical constraint equations for this model are determined uniquely, the 
quantum operators associated with these constraints suffer from the well known ambiguity of the factor ordering, in 
particular in the Hamiltonian constraint. The main purpose of the following will be the motivation of a special choice 
of factor ordering. Technical details of the derivation will be summarized rather briefly. 
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A. Metric representation of the constraint equations 



Let us start with the Einstein Hilbert action for a gravitational field with a cosmological constant A: 

Sbh M = [d 4 x^/=g~ ( 4 ^-2A) (2.1) 



M 



Here M. is the space-time manifold, g = (g^ u ) the 4-mctric on M, g = det(g /MV ) and 4 1Z the curvature scalar of 
the 4-metric. The common prefactor 16 * G of eq. ( |2.l| ) has been avoided by picking units with G — (167r) _1 . As is 

well known, the ADM space-time split and a subsequent Legendre transform with respect to yield the following 
equivalent expression for the Einstein Hilbert action p9| : 



S EH = dt d A x 



^hij-NHc-Ni-H 1 ) , (2.2) 
with H G = G ijk i fr y % kl - Vh ( 3 TZ - 2A) , W = -2 , (2.3) 
where G ijU := -^-= {h ik h jt + h a h jk - h kl ) . (2.4) 



In the transition from ( |2.1[ ) to (2.2) a surface term has been omitted, since it will have no effect on the constraint 
equations. N and iVj are the lapse function and the shift vector, respectively, and the 3-metric (hij) of the spatial 
manifold t — const, is used to raise and lower the spatial indices i, j, fc, . . .. A stroke denotes covariant derivatives 
with respect to the 3-metric, 3 1Z is the curvature scalar of the spatial manifold and h denotes the determinant of the 
3-metric. 

Let us now consider a homogeneous 3-manifold of one of the Bianchi types; then there exists an invariant basis of 
one- forms uj p — U3 p i (x) dx l , such that any homogeneous tensor field on the manifold has spatially constant components 
when expanded in this basis |^,^J. In particular, we have 

du: p = i m™ e qrs u> r A u s , (2.5) 
with a constant structure matrix m pq . In the following we shall be interested in the Bianchi type IX case, where 



the structure matrix is of the simple form m pq = 5 pq . If all the tensor fields occuring in (2.2) are expanded in the 
invariant basis, we arrive at the following expression for the Einstein Hilbert action^: 

S EH = J dt£ EH = J dt J d 3 x u (V 9 hpg - N Ho - N p HpJ , (2.6) 
H = G pqrs n pq n rs - Vh ( 3 K-2A) , H p = 2 e prq m qn h ns n sr , (2.7) 
where h 3 TZ = - (m ps h sp f - m ps h sq m qr h rp . (2.8) 



Here lo = det (w j) in (2.6) contains the only spatial dependence of S EH , which therefore can be integrated out 
explicitcly. Afterwards, the Lagrangian equations with respect to N and N p imply the following set of first class 
constraints 



1 The following expressions are valid not only for the Bianchi IX model with m pq = 5 pq , but even for a general Bianchi A 
model characterized by m vq = m qp . 
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Ho — , H p — 



(2.9) 



with Ho and H p being phase space functions of the fundamental variables h pq and the momenta n pq only.Q If we 
make use of the remarkable identity 



Vh 3 K = ~m pq G mrs m rs 



(2.10) 



which is valid for the general Bianchi type A case, the Hamiltonian constraint Ho may be written in one of the two 
following forms: 



H = [n pq ± i m pq ] G pqrs [n rs qp i m rs ] + 2AVh 



(2.11) 



Let us first restrict ourselves to the kinetic part (which becomes the full Hamiltonian constraint for the Bianchi type 
I model) 



T = n pq G 



pqrs T i 



(2.12) 



describing a particle moving freely on a six-dimensional manifold with coordinates h pq and the (indefinite) supermetric 
Qpqrs ^ w j 1 [ c j 1 jg ^\ ie inverse of G p(?r ,,.[] To quantize such a system, we may employ new coordinates h' pq of the 
minisuperspace such that, at least locally, the supermetric takes a diagonal form, i.e. 



T = tt p<! ' 



t]pqrs ft 



, With Tjpqrs 



-1 



if p = q = r = s = l, 



l(S pr S qs + S ps S qr ) , Otherwise 



(2.13) 



In these 'free falling' superspace coordinates the associated quantum operator is expected to be 



dh 



l"l 



dh' 



(2.14) 



and, transforming back to the /i P9 -coordinates, we arrive at 

VG ohpg 1H dh r , 



(2.15) 



which is easily recognized as the invariant Laplace Beltrami operator on minisuperspace. Here G is the absolute value 
of the determinant of the supermetric G pqrs , which can be shown to be proportional to hT 1 . 

Th ese considerations suggest the following factor ordering for the quantum version of the Hamiltonian constraint 



Ho 



G 



-i h ■ 
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dh 



G Gpqrs 



-i h 
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dh r 



2AVh 



h 2 \fh - 1 G ® 

dhpq y/h Pqrs dh r . 



dhpq \y/h 



T pqrs 



- Vh ( 3 TZ - 2 A) 



(2.16) 



2 According to (2.6) Tr pq are not the canonically conjugate momenta to h pq ; they differ from them by a rescaling factor 
V :— J d 3 xuj. Consequently, the canonical Poisson-brackets read 

{* pq ,h rs } = -±5™ , with S% := i(#<S« + 

and the additional factor y should be carried over to the canonical commutation relations of the quantized theory. There, it 
can be eliminated by a rescaling of Planck's constant Ti. Since we are mainly interested in the quantized theory, we will use 
this freedom to set V — 1 in the following. 

3 More precisely, we require G pqmn G mnrs = 



4 



where we have chosen the signs of im pq in the first line to make the special physical state defined below in ( 2.24 ) 
become an exponentially decaying solution. We end up with a Wheeler DeWitt operator, which, apart from a special 
factor ordering in the kinetic term, contains a quantum correction to the potential, due to the action of the derivative 
operators on the supermetric. Such a term is well-known from the diagonal model, cf. p^| . To put the diffeomorphism 
constraint in more concrete terms, we should be further interested in a /i pg -representation of Tr pq , obeying the canonical 
commutation relations 

[n pq ,h rs ]=-inS Pq . (2.17) 

Moreover, we wish ir pq to be a hermitian operator with respect to the natural auxiliary inner product on the extanded 
Hilbert space, in which the Wheeler DeWitt operator ( 2.1 6| ) is hermitian, 



(*|$) = / d 9 h pq J] KKi - V) y/G <r(M $(M • (2.18) 

J p>q 

Here the components of h pq range over all real values which are consistent with (h pq ) having positive eigenvalues.^] 
As is known from the quantization on a curved manifold, the /i pg -representation of the momentum operator has to 
be corrected with a relative weight factor G 1 / 4 , leading to 

n pq = -i h G~ X ' A - ( J- + J-] G 1 ' 4 =-ih h 1 ' 4 h- 1 / 4 . (2.19) 

2 \Ohpq OhqpJ <J"-(pq) 

Using this /i p9 -representation, we may rewrite Tlo and Ti p in the form 

Ho = (a pq y G pqrs a rs + 2 \Vh = u\ with a pq := h' 1 ' 4 [n pq ~ i m pq ] h 1 ' 4 , (2.20) 



dh 



(sr) 



I L p ■ 



(2.21) 



That the diffeomorphism constraints can indeed be written in the form of (2.21) is checked easily by using the identity 

° h-^ 4 = --h pq , 



dh 



(pi) 



(2.22) 



which implies that there is no additional contribution to 7i p arising from the determinant factor, as one might expect. 
With (2.20 ) and (2.21) we have now nice, self-adjoint operators, which, moreover, form a closed algebra, as we will 
see in eq. (2.3S ) below. 



B. The wormhole state 



As is immediately seen from ( [2.16| ), the wavefunction 



*vrar := ex P 



exp 



(2.23) 



is a solution to the Hamiltonian constraint for A = 0, which moreove r sat isfies a pq fy WH = and therefore solves the 
diffeomorphism constraint (2.21) as well. For the Bianchi IX model, ([2.23) may also be written in the form 



4 The so-defined integration regime for the auxiliary inner product has a non-trivial boundary at h — 0, where at least one of 
the three eigenvalues of (h pq ) vanishes. Consequently, we have to restrict our selves to states ^, which vanish for h — > in a 
suitable manner to assure hermiticity of the differential operators occuring in (2.16). 
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exp 



1 



(Ai + A2 + A3) 



(2.24) 



where X p are the eigenvalues of the 3- metric h pq , and this state is known as the 'wormhole' state from the diagonal 
model.^ The fact that it will occur as a prefactor to all further wavefunctions discussed within this paper suggests to 
perform the following similarity transformation: 



(2.25) 



In this new representation, the transformed operators take the form 

d 



H' = -2ihe prq m qn h„ 



dh 



(sr) 



(2.26) 



n' 



G 



i h ■ 



dh 



2im pq 



i»i 



\/~G Gpn rs i Ti 



dh r 



2k Vh 



(2.27) 



which will be recovered from a very different approach in the next section. 



C. Representation of the constraint equations in Ashtekar variables 

In this section we want to derive the so-called Ashtekar representation of the quantized non-diagonal Bianchi IX 
model using the inverse densitized triad of the 3-metric h pq and complexified canonically conjugate variables thereof. 
The fi rst step in this direction has already been performed in the last section by splitting off the wormhole state 
(2.23). What remains to be done now is to introduce the inverse, densitized triad of the 3-metric h pq , defined via 

h-h™ = e? a -~e\ . (2.28) 



Here and in the following a, 6, c, . . . are flat, internal indices running from one to three; they are raised and lowered with 
the flat metric 5 ab and will therefore always be chosen to be lower indices without any restriction. The introduction of 
a triad has the great advantage that the three metric defined via ( [2.28] ) is automatically positive definite, at least as 
long as the triad is real- valued, and this is favourable for a definition of an inner product on the space of wavefunctions. 
However, as is well-known, we gain three additional degrees of freedom by introducing such a triad, corresponding 
to the three possible rotations in the flat local tangent space. These redundant rotational degrees of freedom are 
accompanied by three additional first class constraints referred to as t he Ga ufi-constraints. 

To construct the generators of these constraints we observe that eq. ( [2.28| ) is invariant under rotations of e p a in the 
local tangent space, generated by 

h d 

J' a = i£abcA pb e p c with A qa :=--—. w ■ (2.29) 

Z U G a 

The J' a satisfy the angular momentum algebra 

[J' a ,J' b }=ine abc J' c (2.30) 



and commute with Ti! and Ti.' . They are the generators of the Gaufi constraints. The general solution of the Gaufi 
constraints J' a ip' = is obviously given by wavefunctions of the form 



5 This name derives from the fact that $ solves the Euclideanized Hamilton Jaco bi eq uati on wi th Euclidean four-geometry at 
large scale parameter. We shall use this name in the following to refer to the states j2.23| ) or (2.24) for A = or its generalization 
for A / 0, without entering a discussion of wormholes, however. 



G 



^'(e p a ) = *'(h pq ) 



(2.31) 



where h pq is a function of the e p a via (2.28). Acting on solutions of the Gaufi constraints the operator d/de p a can be 
written as 



A — — - — ^— = -G g« — 



(2.32) 



where e is the square root of the determinant of e p a . This permits us to rewrite the constraint operators (2.26) and 
( [2.27 ) in terms of the operators e p a and A pa as 



7~tp ■ — 2z £pqr YYl ^sa ? 



(2.33) 



7~^0 " — ^ £abc£pqr ^a^b Q c ' 



(2.34) 



with the operators 



OP — tp 4- — p p — m pq A +- F pqr f u A , A 4- — p p 

a -~ J a i g c a • — <^ x qa ' c abc s^qb -rWc ^ c a 



(2.35) 



The operators Q p a are very convenient and can also be used to reexpress the J' a and 7i' p as 

J a '■= ~£ ^abc A p b Q? c i 



(2.36) 



7~t p — 2i {s pqr e q a Q a A pa J a ) 



(2.37) 



By construction TC and Ti.' coincide with TL' p and Tt' defined in ( 2.26 ) and ( [2.27 ) when acting on the invariant 
subspace spanned by the solutions of the Gaufi constraints (but they extend these operators also to non physical 
states outside this subspace, which are of no interest to us , how ever). The commutation relations of Ti' p and Tt' are 
particularly easy to evaluate in the representation ( |2.33 ), ( 2.34 ) 



[H' p ,H' q ]=ih e pqr m rs H' s , [H' , H' p ] = 



(2.38) 



They imply the same commutator algebra for the 7i' p , Ti.' , and, via the similarity transformation (2.25), also for the 
Tip, Ho 



Summari zing our results so far , we have shown that, assuming that ip'iePg) i s a so lution to the constraints ( 2.29 ), 
( 2.33 ) and ( 2.34j) , the transformed wavefunction ^(/i pg ), connected with ^' via (2.25), is a solution to the constraints 
(2.20) and ( 2.21| ) in the metric representation. In particular, if we are able to solve the more restrictive but simpler 
set of equations 



Q p a i>'(e p a ) = 



(2.39) 



it is clear from the definitions given in (|]||), (|2~36[ ) and ( f2~37|) that we have also found a quantum state of the 
full, non-diagonal Bianchi IX model with a non-vanishing cosmological constant in the metric representation. This 
will only be a special class of solutions, however, because ( ]2.39| ) represents nine conditions, while the Hamiltonian, 
diffcomorphism and Gaufi constraints together constitute seven conditions only. 
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III. COMPARISON BETWEEN THE DIAGONAL AND NON-DIAGONAL BIANCHI IX MODEL 



We shall now adress the interesting question, wether the non-diagonal Bianchi IX model presented in section II and 
the diagonal Bianchi IX model discussed in have the same physical content. On the classical level, it is of course 
unnecessary to distinguish between these two models, because we can use the gauge freedom of the diffcomorphism 
constraints to transform the non-diagonal Bianchi IX model to the diagonal one. When discussing the quantized 
Bianchi IX model, most authors restricted themselves to the diagonal case, i.e. they solved for the diffeomorphism 
constraints on the classical level, and performed the canonical quantization procedure for the effective, 3 dimensional 
system. This approach immediately suggests itself, and one may hope that a quantization of the full Bianchi IX model 
(on a 6-dimensional configuration space) with diffeomorphism constraints imposed on the quantum mechanical level 
should physically lead to the same results. 

In the following, we shall first discuss a very simple example, which immediately shows that this believe is in 
general not true; we will then show in subsection B that the two quantization procedures used for the diagonal and 
non-diagonal Bianchi IX model indeed differ drastically. 



A. An example: The 2 dimensional harmonic oscillator with L = 

Let us consider a well known example, the 2 dimensional harmonic oscillator with unit mass and unit frequency, 
but with the additional constraint that the angular momentum should vanish: 

H = H + NL , 

Wo = g(p?+pl) + 3(g? + g£) ' L = QiP2-q2Pi ■ (3.1) 

Here Ho is the Hamiltonian of the unconstrained system and N is a Lagrangian multiplier. As for the Bianchi type 
IX model, we have the nice property that the constraint is given by a conserved quantity, but L generates gauge 
transformations only if we identify all the directions, in which the effective one dimensional oscillator can move. 

If we firstly quantize this system similarly to the diagonal Bianchi IX case, we have to solve L = on the classical 
level, which is done by 

q 2 = => p 2 = (12 = . (3.2) 
In this gauge we arrive at the effective Hamiltonian 

n eff = \ p2 + ' q = qi ' ( ' 3 ' 3 - ) 
what is easily quantized because it simply describes a one dimensional harmonic oscillator: 

2 2 + °° 
neff=-Y]^ + l q2 > ^l $ >= / d ^*(9)^) ■ (3.4) 

— oo 

Let us now secondly proceed in analogy to the non-diagonal Bianchi IX case. Then we have to quantize first, and 
obtain 

H = H + NL , (V#> = [ d 2 qTp*(q 1 ,q 2 )ct>(q 1 ,q 2 ) (3.5) 

Jr 2 

h 2 ( d 2 d 2 \ I ( d d \ 

with n = -- { M + M ) + - 2 (zl + £) - * = e£-*e£) ■ ^ 

Solving Lip = 0, the wavefunction must be of the form 
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i>(qi,q2) = *(<?) » q = <Jqi + ql 



(3.7) 



However, to ensure that V> is a differentiable function with respect to qi and g 2 , we must require Vt'(q) to be an even 
function in q, 



*(<?) = *(-q) 



(3.8) 



Furthermore, if we compute the action of Ti. on such a solution ^(g), we effectively have 



h 2 1 d d 



1 



t? d 2 



h 2 d 



1 



Heff ~ 2 q dq q dq + 2 q - 2 dq 2 2qdq + 2 q ' 
and the effective scalar product for solutions to the angular momentum constraint becomes 

oo 

(^0) = ^ R2 d 2 gV*( g i ) g 2 )0(gi,? 2 )=27r j dqq 9*(q) $(<?) = . 



(3.9) 



(3.10) 



As a result, we have three important differences between the quantization procedures pointed out above: 



(i) The effective Hamilton operators (3^) and (3J)) differ by a factor ordering term — |^ J^, which becomes singular 



where the coordinate transformation from qi, g 2 to polar coordinates with radius q is not invertiblc. 



(ii) The scalar products (3.4) and ( 3.10 ) contain different weight functions and different integration regimes. 



(hi) For the second quantization procedure, we get a parity requirement (3.8) being absent in the first case, because 
we originally start with differentiable functions on a higher dimensional configuration space. 

All three differences will now be recovered when comparing the diagonal and non-diagonal Bianchi IX model. 



B. Three differences between the quantized diagonal and non-diagonal Bianchi IX model 

To compare our results of section II for the non-d iagon al Bianchi IX model to the diagonal case discussed in jl2) , 
let us try to solve the diffeomorphism constraints (2.21) on the quantum mechanical level. Three special, regular 
solutions to these constraint equations Tip ^(h pq ) — are the invariants of the 3-metric, which read 



T:=Tr(h pq ) , Q := S mn e mpq e nrs h pr h qs , h = det(h pq ) 



(3.11) 



Therefore the general solution to TL p ^ — is any function of these three invariants, 

*(h pg )= X (T,Q,h) , 



(3.12) 



where xiTjQjh) is a differentiable function with respect to its three arguments T,Q,h. We may now express this 
general solution in terms of the three eigenvalues X p of the 3-metric, which, however, are not C 1 (i2 6 )-functions of h pq , 
because of the cubic roots which are needed to express the X p in terms of h pq . Nevertheless, a function ^>(h pq ) solving 
H p W = according to fl3.12| ) actually is C 1 (i2 6 ) with respect to h pq , so 



ip(Xp) = ^(h pq ) = x(Ai + A 2 + A 3 , Ai A 2 + A 2 A 3 + A 3 Ai, Ai A 2 A 3 ) 



(3.13) 



is a regular function in h pq , too. We observe that, as a consequence, tp(Xp) is symmetric under arbitrary permutations 
of the X p , while any wavefunction, which is not symmetric under these permutations will not be differentiable with 
respect to h pq . This symmetry requirement is the analog to the parity requirement (3.8) of our example in subsection 
A. 
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It will be convenient to introduce new variables 

2 



\J XqX r > , with £ p 



(3.14) 



instead of the eigenvalues X p . In the diagonal gauge they play the role of the inverse densitized triad. The inverse 
transformation reads 



An 



fl (Jq <T r 
2 <J n 



(3.15) 



and we arrive at the representation 



*(a p )=^(A p ) = ^(|^: 



(3.16) 



so the wavefunction ^((Jp) is not only invariant under arbitrary permutations of the a p , but in addition invariant 



under reflexions a n 



-cr q ,cr r — > oy. These are necessary symmetry requirements for the wavefunctions 



in the cr p -representation, which were absent in the diagonal case, cf. fll2j. 
The effective Hamiltonian constraint in the <7 p -representation becomes 



— Qi + — 0.2 + — Q 3 

G\ £T 2 03 



*' = , 



(3.17) 



with Q p :— d„d r — d„ 



,d r -a r d q 

1 7 
O q - 



A 

2 a 



v 



(3.18) 



where A := % — , d„ := 



da„ 



(3.19) 



It is possible to show that the very restrictive condition ( [2.39 ) now implies that each operator Q p annihilates ^' 
separately. A comparison with the corresponding result for the diagonal model, cf. eg. ( 2.23) of [12], reveals that, 
apart from a global factor \fo~\ \ we have an additional factor ordering term in ( 3.18| ), which becomes singular 



where two of the a p become identical^} The occurancc of th e sing ularities in these terms and the zeros in the measure 



of ( |3.21 ), see below, reflect the fact that the transformation ( 3.20 ) between the variables h pq in which the quantization 
is performed and the variables {a p ,ip a }, which separate the gauge degrees of freedom from the physical degrees of 
freedom, is not invertible whenever two of the a p take the same value. However, this factor ordering term is just a 
quantum correction, which does not affect the semiclassical limit % — ► 0, and should be expected to appear, cf. (i) of 
subsection A. 

To conclude this discussion, we should remark on the inner product on the Hilbert space of wavefunctions, belonging 
to the two different models. In a first step we want to compute the auxiliary inner product (2.18) for the non-diagonal 
model, which remains for diffeomorphism invariant wavefunctions. Let us therefore rewrite the fundamental variabes 
h pq in the form 



a, 



(3.20) 



where the SO(3)-matrix (£l pr ) depends on three Euler angles tp a . After a subsequent transformation to the a p defined 
in ( 3.14| ) the auxiliary scalar product (2.18) transforms into 



d 3 c 



(<7i er 2 03) 



3/2 



**(0-p)$(t7 p 



(3.21) 



Such a factor ordering term is well-known from earlier work, cf. 122 
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where d£l(ip a ) is the invariant Haar-measure on the group manifold of SO(3) and the variables a p must be positive 
and of a fixed order, e.g. 03 > a 2 > 01 > .[] The Euler angles are pure gauge degrees of freedom and have to be 
ga uge fi xed at some arbitrary value in the physical scalar product, which is simply done by dropping the factor J df2 
in (3.21), which is just a constant prefactor, anyway. 

Remarkably, we get the same result for the effective inner product, if we start from an auxiliary inner product in 
the triad variables e pa , which are defined via 

&pa ' &qa — h p q . (3.22) 



If we impose the simple inner product 



R 



d 9 e P a^*{e pa )4>{e P a) 



(3.23) 



in these triad coordinates we have to eliminate the gauge freedoms with respect to the diffeomorphism and the Gaufi 
constraints. To fix all these gauge freedoms, we may impose the six gauge conditions 



X P a = C P a = 



p^a 



(3.24) 



which require the triad to be diagonal. Then the effective measure in the scalar product (3.23) can be calculated as 
a Fadeev-Popov deter mina nt of the commutators between the generators 7i p , J' a of the gauge transformations and 
the gauge conditions (3.24). If we denote the remaining diagonal elements of e v a by a Pl and, moreover, take into 
consideration all the sy mmetr ies of the wavefunctions with respect to permutations of the a p , we arrive exactly at the 
scalar product given in (3.21). 

In the diagonal Bianchi IX model a natural scalar product is 



(¥|4) 



R 



dad(3+d(3- ¥*(a,/3±)4(a,/3±) 



(3.25) 



because the variables a, /?+ and /3_ introduced in [ fL2[ are the free falling coordinates on minisuperspace for the 
diagonal model. Transforming this to the <7 p -representation, we arrive at 



(¥4) oc 



cr„>0 



d<7i d<j2 das 

<J\ (72 03 



¥> P )4K) , 



(3.26) 



where the three <7 p -integrals have to be performed along the positive real axes now. Obviously, the additional factors 



(7p — a occuring in (3.21) are absent in the diagonal case, i.e. the measures in the scalar product of the diagonal and 



the non-diagonal model differ, as was already expected when discussing point (ii) in subsection A. 

We conclude that quantization and gauge- fixing are essentially not interchangable for the model under investigation. 
The quantization procedure followed in section II seems preferable because it seems much nearer to a quantization 
of the full, inhomogeneous field than the methods used for the diagonal case. Moreover, all four constraints TC P are 
quantized and in this sense treated in a similar way, as one would expect from the viewpoint of general covariance. 

To complete our discussion of the inner product on the sp ace o f wa vefun ctions , we should give the physical inner 
product, which is obtained from the auxiliary inner product (3.21) or (3.26) by eliminating the last gauge degree of 
freedom corresponding to the Hamiltonian constraint H.q. 



C. The physical inner product on the Hilbert space of wavefunctions 

To shorten the present paper, we want to postpone a detailed motivation and derivation of the physical inner 
product to a future paper 118[. There, we will also investigate the normalizability of the wavefunctions constructed in 



7 The transformation (3.20) turns out to be invertible only for a fixed order of the X p (or, eqivalently, of the a p ) 
permutations of the A p may be realized by suitable rotations with SO(3)-matrices Q pr . 



because 
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section IV, and it will turn out that the physical inner product uniquely determines one normalizable state in both of 
the two cases A = and A > 0. In this subsection we only want to give the final form of the physical inner product 
on the space of wavefunctions and some main ideas how to arrive there [p3|~p5| . 

The effective inner product for a constrained system with constraint T — is generally obtained by firstly choosing 
a suitable gauge condition and secondly calculating the following expectation value with respect to the auxiliary 
(or 'kinematic') inner product: 

(¥||*W = <¥||J|a(x)|*>««™ , with J:=i[T, X ] ■ (3.27) 



Here J is the so-called Faddev-Popov determinant. This definition makes sense whe never T is a first order derivative 



operator on configuration space, because then \ J\ will be a positive real number and ( 3.27 ) defines a positive hermitian 
product. However, in the present case T = H.o is a second order derivative operator, and we have to deal with two 
problems: 

(i) At first sight, |J| seems to be a mysterious quantity, because J itself is a first order derivative operator on 
configuration space. However, since J is a self-adjoint operator, we may decompose any state \t in terms of 
eigenstates of J to arrive at a very natural notion of |J|. This operator |J| then is itself self-adjoint, and, 
moreover, positive. 

(ii) In general, |J| and 5(x) will not commute, so \<&} p hys will neither be hermitian, nor positive. This would 
give a rather miserable physical inner product. Surprisingly, there exists a small set of gauge conditons \ such 
that the gauge condition commutes with the Faddev-Popov determinant, i.e. such that x solves 

[[T,X],X] =0 • (3.28) 



It seems that this set of gauge conditions does not only exist for the case under study, T = Ho, but for any 
physically interesting case of constraints which are described by second order derivative operators. Using these 
special gauge conditions the physical inner product (3.27) obviously is poitive and hermitian. 



Having solved these two problems, we are now ready to perform the last step in the construction of the physical 
inner product for the non-diagonal Bianchi IX model. It is easily seen that the gauge condition 



X(<7l, 02,03) 



2, 3 

-In 

n 03,0 



(3.29) 



is a solution to eq. ( |3.28 ). It fixes the variable 03 to have a specific value 03 
transformations the associated inner product finally takes the following form: 



03.0- Performing suitable coordinate 



d/3_ / dv V* 



\ lh ir $ 

1 ov 1 



with 



(3.30) 



(af - a\) (0. 



05 



01 02 03 



1/2 



v = a 



3/2 

- ) V "! a 2 CT 3 



(3- = In— . (3.31) 
2^3 01 



This inner product has now an intrinsic 'time'-parameter 03, which must be given to compute probability distributions 
with respect to /3_ and v, the remaining physical degrees of freedom. 



IV. SOLUTIONS GENERATED FROM THE CHERN-SIMONS FUNCTIONAL 



A. Ashtekar representation of the Chern-Simons state 



As pointed out in section II, we can find a special exact quantum state for the non-diagonal Bianchi IX model with 
a non vanishing cosmological constant, if we are able to determine a function ip' lying in the kernel of Q p a defined 
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in (2.35). To proceed in this way, let us consider the Fourier transformed version of (2.39), where it is essential to 
perform not a customary, but a generalized Fourier transform introduced in [fl2|| : 



<//(e p J oc / d 9 A pa exp 

J£9 



A pP 

■"■pa c a 



*(A pQ ) 



(4.1) 



Here E 9 C C 9 is any 9 dimensional manifold in the 18-dimensional space of the complex A pa , which allows for partial 
integration without getting any boundary terms. Then ip'(e p a ) is a solution to ( |2.39 ), if ^(A pa ) solves the following 
set of equations: 



i \ 8 

■,p<i A p pqr f t. A h A A — 

1 -"qa ~ £ tabc f±qb "-rc ' c\ a 

Z 2 OAp a 



®(A pa ) = 



(4.2) 



A special solution to (4.2) is the Bianchi IX restriction of the well known Chern-Simons functional: 



* cs = exp 



i ( - A pa m pq A qa + 2 det(A pa )) 



(4.3) 



In the following we shall be intere sted in the transformation of this state to the metric representation, i.e. in the 
evaluation of integrals of the form (|4.l| ). Topologically different choices of E 9 (i.e. choices which cannot be deformed 
into each other without passing a singularity) will lead to different states in the metric representation. 



B. Metric representation of the Chern-Simons state 

It will be convenient to introduce new variables 

K V ■= ^ 7.V — ^1 



(4.4) 



Then the wavefunction in the metric representation takes the form 



*' cs cx / d 9 A pa exp 
Js 9 



^ ( K p a A P a — 2 Apa rn pq A qa + det(A pQ ) 



(4.5) 



where suitable integration manifolds E 9 remain to be determined. Here and in the following, constant prefactors, and 
also prefactors depending on A, are absorbed in the proportionality sign 'oc', since such a factor will always remain as 
a freedom of the wavefunction. Surprisingly, six of the nine A pa -integrals turn out to be Gaussian integrals and can 
be evaluated analytically, namely the integrals with respect to An, A\ 2 , A13, A21, A22 and A23. Lengthy calculations 
finally give the result 



dx dy dz 



1 



exp 



x{- r +2K3 ' r + — 



\2 



(4.6) 



where we have introduced the abbreviations r — (x,y,z) :— (A 31 , A 32 , A 33 ) and n p :— (k p 1; k p 2 , k p 3 ). It is clear that 
the 3-direction is distinguished by the order in which the integrals are performed. As we know from section II, Sb' cs 
according to (4.6) is automatically Gaufi- and diffeomorphism invariant, because it is a solution to (2.39). This means 
that ty' cs can only depend on the eigenvalues A p of the three- metric h pq , or, equivalcntly, on the a p introduced in 
( |3.14 ). So we are free to choose a diagonal gauge for the triad in (4.6), i.e. we may take 



(« P a) = diag(«i, k 2 , k 3 ) 



With K3 > K 2 > Kl 



(4.7) 



where k v 



A 
12 



\J Aq A r 



(4.8) 



This special ordering of the three n p (or, equivalcntly, of the A p ) is always possible, because permutations of the 



13 



diagonal elements of a diagonal metric h pq correspond to diffeomorphisms, which leave the wavefunctions invariant. 
Moreover, such an ordering is even necessary to have a unique mapping (k p q ) — > (ki,K2, K3) under GauB- and 
diffeomorphism transformations, and we will adopt the convention (4.7) throughout the following. In the gauge (4.7) 
the integral (4.6) takes the form 



dx dy dz 

oc / — 2~exp 

/ S 3 1 - r 2 



2 K3 z - 



2 K\H2 z - 



2 2 

«2 y 



(4.9) 



where r 2 := x 2 + y 2 + z 2 . We will see in sections V and VI that there are choices for £ 3 in (4.9) leading to 
asymmetric wavefunctions under (formal) permutations of the k p ^ However, as we know from section III (cf. after 
eq. ( 3.16| ) with k p — ^ o~ p ), only wavefunctions which are symmetric under permutations of the k p and under reflexions 

of interest to us, and this will be an important restriction to select the physically 

interesting wavefunctions. 



V. ASYMPTOTIC FORMS OF THE CHERN-SIMONS INTEGRAL 



In order to get information about possible integration manifolds £ 3 which can be used in eq. ( |4.9| ) we begin by 
discussing the asymptotic behavior of the Chern-Simons integral in several asymptotic regimes. It will thereby become 
obvious that we deal with five linearly independent solutions, and the integration contours will be given in section VI 
afterwards. 



A. The semiclassical limit h — > 



Let us first of all examine the semiclassical behavior of the solutions described by (|4.9|) . Surely, the semiclassical 
limit could have also been discuss ed b y starting from the nine dimensional integral (|4.5| ), but, of course, an expansion 
of the three dimensional integral (|4.9| ) is much simpler. The saddle point form of d4.9|) is displayed nicely by writing 
it in the form 



dx dy dz 

1_„2 eX P 



3F 



(5.1) 



with F 



+ 2k 3 z ■ 



2 K1K2 z — k\ x 2 



2 2 

4 y 



1 



(5.2) 



where we have inserted A according to ( 3.19| ). Since the n p defined in ( |4.8| ) do not depend on h, we approach a 
Gaussian integral in the limit Ti — > 0, A fixed, which has to be evaluated at one of the saddle points of the exponent. 
However, which of the possible saddle point contributions arises for the integral under consideration is determined by 
the integration surface S 3 and requires a detailed discussion of the contours of steepest descent. Here we only want to 
give the possible asymptotic results for h — > which may be realized by suitable choices of the integration contours. 
The integration contours are discussed in section VI. 

The saddle points of the exponent are obtained by solving the equations 



dF 

dx 



= 



dF 
dy 



= 



dF 



= 



(5.3) 



One can show th at f or <7i,er2,<73 pairwise different, i.e. 
only solution to (|]^) is given by 



in particular on the sector 173 > 02 > <J\ > of interest, the 



x = y = , (z - k 3 ) (1 - z 2 ) 2 = (ki + k 2 z) (k 2 + k\z) . 



(5.4) 



In complete analogy there are eigenstates of the effective Hamiltonian (3.9) of a 2-dimensional harmonic oscillator with 
L = 0, which do not obey 'i'(q) = g), but these are of no physical interest. 
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A comparison with the diagonal model p2[ reveals that for x — y = the exponents of ( |5.l| ) and of the one dimensional 
integral (4.17) of JL2| become identical. Consequently, also the saddle point equations with respect to z coincide, cf. 
eqs. (5.4) and (4.7) of Jl2|. As a result, the semiclassical actions of the wavefunctions approached in the limit h — > 
must be the same, since they are given as the saddle point values of the integrand's exponents. This is of course 
exactly what we expected to happen, because the Wheeler-DeWitt equations for the classically equivalent diagonal 
and non-diagonal Bianchi IX model give rise to the same Hamilton Jacobi equation in the limit ft — > 0. However, the 
quantum corrections to this leading order behavior are rather different, firstly because of the different dimensionality 
of the integrals, and secondly because of a new prefactor to the exponential function. Performing a Gaussian saddle 
point approximation, we get for the non-diagonal case 



, TWO 



(-2^A) 3 

dr a dr b 



det 



F/\ 



d 2 F d 2 F d 2 F 
dx 2 dy 2 dz 2 



e F ' x 
l-r 2 



(5.5) 



Here the second proportionality stems from a special property of the exponent F, namely the fact that the mixed 
derivatives of F vanish for x = y = 0. Comparing this with the asymptotic results obtained in the diagonal case (L2J, 
we find that corresponding results differ only by a factor 

_ j_ 

(5.6) 

x=y—0 : z—z s 



7 := 



, x d 2 F d 2 F 



dx 2 dy 2 



Since we already know the asymptotic results for the diagonal Bianchi IX model, we only have to expand this additional 
factor at the saddle points to get the quantum corrections for the non-diagonal Bianchi IX model. As in the diagonal 
case, we therefore have to deal with a saddle point equation of fifth order, so analytical results are only available in 
additional asymptotic regimes. Due to the existence of five saddle points one has five linearly independent solutions. 



1. The limit A -> 



One asymptotic regime which allows for an expansion of the five saddle points is th at of small cosmological constant 
A — > 0. Since we are in the additional limit h — ► 0, the quantity A = defined in ( 3.1 9| ) tends to zero, too. In this 
limit one saddle point is given by 



A^O 



1 ^ x ' 



(5.7) 



The contribution from this saddle point yields 



W^g tx const. , (5. 



and therefore approaches the well-known wormhole statej^j 
Four further saddle points asymptotically lie at^| 



A^O A , , A^O A , , 

z - -l±-(a 2 -ai) , z - l±-(<7 2 +<7i) , (5.9) 



giving rise to an asymptotic behavior 



9 We r emind the reader that the prime at ^' C s denotes the fact that this wavefunction has to be mutliplied with ^> W h defined 
in (2.24) to become a solution of the Wheeler DeWitt equation in the metric representation. 

10 The following expansions are only valid for 02 > 01, but we restricted ourselves to this case anyway. 
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A^O 

, cx 



-C3T( a '2-<Tl) 



A^O 

, CX 



+ <7 2 )(ct 3 T o-i)(ct3 T cr 2 ) 



(5.10) 



respectively. These solutions are the generalizations of the diagonal analogues (2.25) and (2.26) given in |12|. The 
singularities which oc cur in the denominators of the asymmetric solutions of (5.10) are cancelled by the weight function 
of the inner product ( 3.21 ) and therefore do not constitute any physical problems. 



2. The limit k — > oo 



A second regime which allows for an analytical expansion of the semiclassical limit h — » is the case n — > oo, where 
the 'mean' k is understood as K := (k% ^2 K3) 1 / 3 . The three k p where defined in eq. (4.£). For a fixed cosmological 
constant, this is the limit of large overall scale paramter. In this limit, two of the saddle points behave according to 



± 



K3 



(5.11) 



and turn out to be complex-valued. This gives rise to a complex action of the wavefunction, describing a Universe 
with a Lorentzian signature of the 4-metric, cf. pjfl. The leading order behavior for these solutions is given by 



exp 



±- 



4i AC 3 / 2 

A 



(5.12) 



for a more detailed asymptotic result we refer to |12| . Three further saddle points have an asymptotic expansion 



K — ► 00 K\ 



k - — * 00 K>2 



1 + k 3 , 



K2 Kl 

and lead to wavefunctions of the form 

*' cs K oT [(4-4) (4-^)T^ 4/x , *' c /cT [(4-4) («g -«?)]"* e«?/* 



(5.13) 



*' c /*°° [(4-4) (4-4)]'^^ 



(5.14) 



respectively. These saddle point contributions are known from the diagonal model for the 'asymmetric' solutions 
introduced there; however, here we have additional prefactors, which become divergent on the symmetry lines a p — 
Vqi P 7^ <7j but again these singular terms are cancelled by the measure of the scalar product (3.21). 



B. Solutions to the vacuum model approached in the limit A 







In order to get solutions of the non-diagonalnon-diagonal Bianchi IX model for A — > (without taking the semi- 
calssical limit Ti — > first) we now want to discuss the behavior of the Chern-Simons integral in this limit. While one 
of the solutions approached for A — > again turns out to be the wormhole state (2.24), the other four solutions are not 



given by the asymptotic results ( 5.10 ), as one might think at the first sight. The reason for this is the fact that we are 
not allowed to perform the usual Gaussian saddle point expansion, since the exponent, and in particular the prefactor 
to the exponential function, become singular in the limit A — > 0, where r approaches ±1. Comparing the situation 
with the diagonal model, where a saddle point expansion in the same limit actually was allowed, one might ask for 
the difference between the two models that forbids for an analogous expansion in the present case. The answer to 
this question lies hidden in the prefactor to the exponential function: while we dealed with an integrable square-root 
singularity in the diagonal case, we here have to integrate into a singularity of first order, even in the case A = 0, and 
this requires much caution. In fact, it turns out that the calculation of the limit A — > of the Chern-Simons integral 
is very subtle and worth to be discussed in a separate section in appendix A. Moreover, we show in appendix B that 
the solutions derived in appendix A can be written in another, very nice and compact form, which reads 
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r 3 

/ dse s ]J( S -s u )-i 



(5.15) 



Here the integration contour C e is one of the four curves shown in figure [l| and the s v are special sums of the a p , 
defined by 



«o = o"i + c 2 + cr 3 , si = cti - a 2 - 0-3 , s 2 — -cri + a 2 - a 3 , s 3 = -o\ - a 2 + cr 3 



(5.16) 



For a detailed discussion of these wavefunctions and a comment on their symmetries with respect to permutations of 
the a p we refer to appendix A and B. 

Altogether, wc find in the limit A — * again five linearly independent solutions associated with the Chcrn-Simons 
wavefunction (4.5). 



C 



12 



A 



-*n¥t 



c 



123 



WWW 



c. 



0123 



FIG. 1. Integration curves in the complex s-plane which lead to linearly independent vacuum solutions. A dashed line 
indicates that one has to evaluate the integrand in the second Riemannian branch after having crossed one of the cuts. The 
two cuts are represented by wavy lines. 



VI. THREE DIMENSIONAL INTEGRATION MANIFOLDS FOR FIVE EXACT SOLUTIONS 

In this section we want to define explicitely five integration cont ours S 3 for the integral mM), leading to five 



linearly independent solutions of the Wheeler DeWitt equation ( 3.17 ). To have a better view ove r the symmetries of 
the solutions, let us, for a moment, go back to the nine dimensional integral representation ( |4.5| ). It can be verified 
easily that the two different integration surfaces 

S 9 ± := {(4 pa )eC 9 |A pa eRe ±i7r / 6 } (6.1) 



lead to a finite integral, because the cubic terms in the exponent of (4.5) remain purely imaginary, while all the 
quadratic terms have a ne gative real part. Moreover, since the only border of this surfaces lies at infinity, it is 
clear that the integral ( |4.5| ), performed over Sj_, must lead to a solution of the non-diagonal Bianchi IX model. In 
particular, we then know that these wavefunctions are Gaufi- and diffeomorphism invariant, and we may again choose 
the diagonal gauge (4.7) for the matrix (n p a )- Then (p~5|) takes the form 



cs,± 



d 9 A 



exp 



Ki An + k 2 A 22 + k 3 A 33 - - A pa m pq A qa + det^,,) 



(6.2) 
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If we now consider a formal permutation n\ <-> k 2 in (6.2), we can re-establish the original integral by a suitable 
transformation in the ^4 pa -space, namely An <-> A 22 , ^32 «-* ^.13, ^23 ^.31- Since all integration variables A pa 
are integrated along the same axes in the complex A pa -pla nes, the integration contour remai ns u nchanged under this 
coordinate transformation and we regain the integral (6.2). Furthermore, the two integrals ( |6.2| ) are invariant under 
a substitution K\ — ► —K\, k 2 — ► —k 2 , which may be seen by a transformation A\ a — ► — A\ a , A 2a — y — A 2a in the 
A pa -sp&ce. Thus we conclude that both wavefunctions are completely symmetric under arbitrary permuations of the 



k p , and also under reflexions n p 



a symmetry property which is posessed only by two of 



the five solutions which are generated by the Chern-Simons f unct ional 



Performing the six Gaussian integrations which lead from (4.5) to (4.9), we reach the following three dimensional 
integration surface for the integral (4.9): 



S 3 ± 



(6.3) 



In the following it will be convenient to introduce the two superpositions 



vff' . oc vEr' , 

* WH.A * cs,+ 



OS,- 



(6.4) 



and also new integration variables 



x — r yl — T] 2 



cos ip , y 



= r yl — if- 



sin (p , z = rr) 



(6.5) 



for the x, y, z-integral, which may be understood as spherical coordinates with rj = cos-d. For the transformation (6.5) 
the volume clement simply transforms like dx dy dz = r 2 dr drj dtp, and the condition r 2 = x 2 + y 2 + z 2 is obeyed even 
for com plex values of r, 77, ip. With this new parameterization, the integration surfaces for the two solutions defined 
in (6.4) can now be written in the following form: 



= {re C 3 I - 1 < rj < 1, < ip < 2ir, r e R + i e) 



e > 



(6.6) 



s oi23 = {reC 3 I - 1< rj < 1, < ip < 2tt, reK(l,R)} , R < 2 



(6.7) 



In ( |6.7D K(1,R) denotes a pole integr al a round r = 1 with radius R, cf. appendix A. Let us try to interpret these 
integration surfaces: in the first case (3J3) we simply integrate over a sphere with radius r in the complex r-space; 
however, this radius is not chosen real- valued as usual, but slightly imaginary to avoid the singularities of the integrand 
lying at r — ±1. The integration manifold Sq 12 3 also describes an integral over a full sphere, but now the radius 
des crib es itse lf a circle round the singularity at r = +1 in the complex r-plane. While the two parameterizations given 
in (3J3) and (6/7) will prove useful for an asymptotic expansion in the limit A — > 0, which is performed in appendix 
A, they now allow for a further study of the scmiclassical behavior of the corresponding solutions. 

If we consider the wavefunction v E , o 12 3 a m ^ ne semiclassical limit Ti — > 0, we get the dominant contributions to the 
integral for rj = ±1, so that the integrand becomes independent of ip. The remaining r-variable then coincides with 
the invariable, which was the only integration variable in the diagonal Bianchi IX model discussed in p^] . For a 
discussion of the curves of steepest descent in the z-plane we refer to this earlier work. However, if one has finally 
deformed a desired integration contour into the curves of steepest descent, there remains one remarkable difference 
between the two models: while in the diagonal case we had a cut in the complex z-plane for |Rez| > 1, this cut is 
absent in the present non-diagonal case, owing to the different prefactors to the exponential function. So integrals 
along parts of the real axis of the z-plane, which cancelled in the diagonal case, may now be different from zero, and 
vice versa. A detailed discussion of the curves of steepest descent and the result ing a symptotic contributions finally 
gives the following result for the wavefunctions 4", being connected with via ( 2.25 ): 

In the semiclassical limit Ti — > the wavefunction ^0123^ shows a no-boundary behavior for a — * 0, i.e. the 4- 
manifolds which correspond to the semiclassical trajectories in minisuperspace are regular at the point a = 0. For 
a — > 00 , the integration surface Sq 123 picks up several saddle point contribution s: th e two complex conjugate saddle 
points (5.11) contribute as well as the two real- valued, negative saddle points in (5.13). These latter two saddle points 
are the reason why v I'oi23 A is not normalizable in the physical inner product (3.30), as will be shown in a future paper 

PL " 
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The solution 'EVh.A is a generalization of the wormhole state (2.23) for a positive cosmological constant A > 0. For 



a — > oo, we again gather contributions from the two c omple x conjugate saddle points ( 5.11 ), and now, in addition, 



from the real- valued positive saddle point z ~ 1 + K3 of (5.13). In |18J we will see that these latter three saddle points 



give rise to finite co ntrib utions in the physical norm of the wavefunction. We will show this by evaluating the physical 
norm according to (3.30) in the semiclassical limit h — > 0. As a result, ^wh .a tu rns out to be the only Chern-Simons- 
like quantum state, which is normalizable in the physical inner product ( |3.30 ) and symmetric under permutations 



and reflexions of the a v . The existence of a normalizable Chern-Simons-like state is in remarkable contrast to earlier 
work on special spatially homogeneous models discussed by Marugan [^6| . 

Finally, we still have to define three further integration contours Y? to obtain the complete set of Chern-Simons like 
solutions. The construction of the vacuum solutions in appendix A suggests three new possibilities to create linearly 
independent integrals. One of these further integration surfaces can simply be written in the form 

Y? 12 = {feC 3 I -oo< r) < -1, < Lp < 2tt, reK(l,R)} , (6.8) 



while the other two have a very complicated form and may be read off from the representations (A2C) and (A21) in 
appendix A. That these latter two integration manifolds, despite of their strange ?y-junctions at — ci/cra and — (72/ci, 
give rise to analytical solutions of our model even for A =/= is shown in appendix C. However, these solutions do not 
obey the necessary symmetry properties which were derived in section III, and are therefore of no further interest to 
us. 



One may now finally ask, how the four solutions defined in ( AlS ) to ( A21 ) of appendix A are connected to the 
'asymmetric' solutions of the diagonal model discussed in |12fl. An investigation of the semiclassical behavior, 
together with a discussion of the singularities on the lines a v = a q , p ^ q, shows that the solutions *&i2i ^123 an d 
^0123 are actually the sums of these states as indicated through the choice of their indices, i.e. we have 

f ' 12 = + * 2 , *' 123 = + f 2 + *' 3 , * 0123 = % + % + * 2 + * 3 ; ^3 > 02 > O-l . (6.9) 

By inverting these relations with respect to one may further define asymmetric solutions corresponding to those 
of the diagonal model, with the same symmetry properties that have been pointed out in fl^| . 

VII. DISCUSSION AND CONCLUSION 

In this paper we have examined the transformation connecting the representations of quantum general relativity in 
metric variables and in Ashtekar's variables for the special case of spatially homogeneous but anisotropic space-times 
of Bianchi type IX with non-diagonal metric tensor. While classically the non-diagonal case and the diagonal case 
are equivalent (in the absence of matter) due to the freedom of gauge-fixing, there is a subtle difference quantum 
mechanically, because the steps of gauge-fixing and quantization, in general, do not commute. This was explained via a 
simple example in section III A. The example also made clear that gauge- fixing after quantization is preferable because 
all symmetries are then implemented automatically and democratically. Once the two steps have been completed in 
this order it is then also clear a posteriori how to proceed correctly in the quantization after the gauge-fixing has 
been performed on the classical level. Performed in this way quantization and gauge-fixing of course commute by 
construction. 

For the non-diagonal Bianchi type IX model we have considered in detail the transformation of the Chern-Simons 
state from Ashtekar variables to metric variables. The Chern-Simons state, including its limit for vanishing cosmolog- 
ical constant, deserves a thorough study because it is undoubtedly the most important exact solution of all constraints 
of quantum general relativity found up until now. This is so because, unlike all other exact solutions, it describes a 
well-defined non-degenerate space-time in its classical limit and also because it makes an obvious connection between 
quantum general relativity and topological field theory |t],|8| . Indeed it is remarkable that the general Chern-Simons 
state in Ashtekar variables makes no reference to metric concepts on the space-like 3-manifold. It should be noted 
however, that writing down a physical state in Ashtekar variables does not yet define it completely, because 'reality 
conditions' must still be imposed before a physical interpretation can be attempted. Imposing the reality conditions 
is a very nontrivial task and, as our results indicate, might not have a unique solution. We circumvent this problem 
completely by transforming back to the real metric representation before applying a gauge-fixing and giving a physical 
interpretation. 

The Ashtekar variables and the densitized inverse triad form canonically conjugate pairs. Hence, for our Bianchi type 
IX model the generalized multidimensional Fourier transformation we discussed in sections IV, V, VI can be used to 
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transform the Chern-Simons state to the metric representation. This is a generalized Fourier transformation because 
neither the integration contours, which are here 9-dimensional manifolds, nor their boundaries are fixed a priori, 
except for the condition that partial integration with the Chern-Simons state under the integral must be permitted 
without contribution from the boundaries. The boundaries are therefore determined entirely by the singularities of 
the Chern-Simons state. On the other hand, for fixed boundaries the different integration manifolds one can find may 
be deformed without changing the result, or may be combined by first running through one integration manifold and 
then through other inequivalent ones, leading to linear combinations of the physical states defined by each integration 
manifold separately. 



In eq. (4.9), with the five integration manifolds £ 3 given by eqs. (6.6), (3.7), (6.8), (A.2C), (A.21), we have obtained 
integral representations of five exact solutions to all constraints of quantum general relativity and we studied various 
limits of these solutions in section V. In the leading semiclassical order the result for the non-diagonal and the diagonal 
model is the same, as one would expect from the classical equivalence of both cases. However, differences appear 
already in the next to leading order, obtained in the semiclassical expansion of our results in section V A. It should 
be noted that even studying the next to leading semiclassical order of a physical state in quantum general relativity 
is rather nontrivial. It requires to address operator ordering ambiguities in the Hamiltonian, which we have done 
in section II A, to take proper account of quantum corrections in the Hamiltonian implied by symmetries, which 
we have done in section III B, and to apply a complete set of gauge-fixing conditions including a fixing of the time 
reparametrization symmetry, which we have done in section III C. Only then is the next to leading semiclassical order 
of a physical state unambiguously defined. 

Only two of the five solutions we constructed satisfy the complete permutation symmetry between the three main 
axes of the Bianchi IX 3-geometry, which, as shown in section III B, is implied by the quantization of the non- 



diagonal model. These are the states corresponding to the integration contours d6.6| ) and (6.7). Their semiclassical 
limits identify them as a wormhole state and a Hartle-Hawking 'no-boundary' state, respectively. It is remarkable 
that the same Chern-Simons state in Ashtekar variables can yield such diverse states depending on the choice of the 
integration manifold. This is a striking example that a physical state in Ashtekar variables is not yet defined before 
further conditions (here the choice of integration contours) are specified. 

In the Hartle-Hawking proposal for the semiclassical initial condition of the classical evolution of the Universe it is 
not required that the 'no-boundary' state is a normalizable vector in a Hilbert space. Rather the condition by which it 
is defined at least semiclassically is that the Euclidean 4-geometries defined by the semiclassical wavefunction filling in 
the 3-geometries on a space-like slice are regular for a — ► 0. A further requirement for the Hartle-Hawking state is to 
give a well-defined probability distribution at the semiclassical caustic surface where the semiclassical evolution, given 
by the wave function, switches from a Riemannian ('Euclidean') to a Pseudo-Riemannian ('Lorentzian') space-time. 
This requirement is met by the Hartle-Hawking 'no-boundary' state obtained here, and only by this state among the 
five Chern-Simons-like solutions, as was shown already in Jl^ |. 

The wormhole state, on the other hand, one expects to be a normalizable vector in a Hilbert space with a well- 
defined scalar product, in which other state vectors describe e.g. excited states of the wormhole. Without going into 
details, which will be given elsewhere |Q , we have defined such a scalar product in section III C by gauge-fixing all 
gauge symmetries. It turns out that the wormhole state is, in fact, normalizable in this scalar product, while the 
'no-boundary state' is not. Thus imposing either the 'no-boundary' condition or normalizability in the Hilbert space 
of physical states equipped with a fully gauge- fixed scalar product, one finds in each case a unique but mathematically 
and physically vastly different state as a metric representation of the Chern-Simons state. 

One may hope that the non-diagonal quantization procedure presented here might be applicable even to the full 
inhomogeneous case, because all constraints are treated quantum mechanically in a similar way, and a solution to 
the quantum constraints is available in terms of the general Chern-Simons functional. Work in this direction is in 
progress. 
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APPENDIX A: THE CHERN-SIMONS INTEGRAL IN THE LIMIT A -> 



In the following we want to examine the possible limits A — > of the Chern-Simons integral given in(4.9). First of 
all, we shall briefly recover the wormhole state, which is approached for the choice £ = £^- H given in 7|j). Using 
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this integration manifold, we may write 



+ 1 27T +oo+ie ^ 

« J dr) J dip J ex P 



(Al) 



— oc+ie 



where we made use of the spherical coordinates in the (x, y, z)-space, which were introduced in (3.E). Easy estimates 
show that for A — > just an infinitesimal region around r — contributes to the inner r-integral, so we have the 
following asymptotic behavior:^ 



+ 1 27T +e +oo 

% VH *«° «P 47r A 3 / 2 I d^ 2 e"« 2 = 2 (.A) 3 / 2 



const. 



Vac . (A2) 



According to the transformation rule (2.25) we obviously approach the wormhole state ^ WH in the limit A — > 0, which 
was defined in (gjg). 

Let us now turn to the much more complicated cases, for which the 'wormhole saddle point' x = y = z = 0is not 
passed through by the integration surface S 3 . For such solutions, we shall try to set not all borders of S 3 at infinity; 
instead we shall make use of the existence of two further singularities of the integrand with respect to r, namely 
those at r = ±1. If we integrate into these singularities in a suitable manner, we expect to create further solutions of 



our Wheeler-DeWitt equation (3.17), because boundary terms generated by partial integrations will vanish at these 
borders. In the following it will be in fact sufficient to consider only one of the two singularities, say r = +1, since 
any r-integral in the neighbourhood of r = — 1 can be mapped onto a region around r = + 1 by a transform r\ — > — r\ 
in the ^-integral. This is due to the fact that the r/-,r-dependence of the integrand in (4.9) is given in terms of r 2 ,ri] 
and rj 2 only, cf. (A4) below. Thus we will be interested in integration countours for the r-integral, which have one 
end point at r = 1. It turns out that, for a positive cosmological constant A, there are only two curves of interest: 
firstly, we may perform our r-integral along the real axis from r = 1 to r = +oo; secondly, we can also consider a 
pole integral, encircling the singularity at r = 1 in the mathematically positive sense. To shorten our notation, we 
should establish the following nomenclature: a circle with radius R, which is centered at r and passed through in the 
mathematically positive sense will be denoted by K (tq, R). Then we can write the integrals of interest in the form 



r <■ 



dry / dip 



1 dr 



1 -r 2 



exp 



with 



(A3) 



G 



2n 3 ri]+ -(ft 2 



1 (k 2 + k 2 )(1 + r2 V 2 ) + 4 «i K2 r n + r 2 (k 2 — k\) {n 2 — 1) cos2ip 

2 1-r 2 



(A4) 



where C is one of the two contours ]l,+oo[ or K(1,R), < R < 2. The equality sign in (A3) indicates that, from 
now on, all prefactors of the integral will be taken into account. To proceed in our calculation of the limit A — > 0, we 
can uniformly estimate those parts of the integral ( |A3| ) , which lie outside an e- neighbourhood of r = 1 by a function, 
which vanishes exponentially for A — > 0. Therefore, just an e- neighbourhood of r = 1 contributes to the integral QA3|), 
where e may be chosen arbitrarily small. As a consequence, all terms of the integrand, which remain regular at r = 1, 
can be substituted by their value taken at r = 1, and only the r- dependence of singular parts remains to be taken 
into consideration. Thus we arrive at once at 



dip 



dr 



exp 



X (l-r) + 7 



A a 2 — 2 ab cos 2ip + b 2 



1 



(A5) 



where we made use of the definitions 



Cl + 0- 2 , . <7i - <7 2 , , 

a := (1+r?) , 6:= (1 - rj) 



(A6) 



As h is kept fixed, also A = cf. (3.19), tends to zero in the limit A — » 
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Here C £ is either equal to ]1, 1 + e[ or given by K(l, e). After the transformation £ = 4(r — 1) we arrive at 

£ a 2 — 2 ab cos 2ip + b 2 ~ 



2i 



(A7) 



where now C is the positive real axis or K(0,R). The two possible ^-integrals can now easily be performed, leading 
to 



A ~° / d77 e CT3 " / d(^ 




^/a 2 - 2 a6 cos 2^ + 6 2 ,C' = K(0,R) 



y/a 2 - 2 ab cos2p + b 2 ,C'=]0,+oo[ 



(A8) 



Up to now we have taken into acount all freedom in the choice of the integration cont ours ; any other choice fo r th e 
r-integration contour would have given a linear combination of the two results given in ( A8) or a constant as in (|A2| ). 
While the ^-dependence of flA8|), hidden in a and b, is to complicated to be integrated out analytically, there are 
several possibilities to evaluate the ^-integral, which become more transparent in the new variable u = cos 2ip. With 
the abbreviation X := a 2 — 2abu + b 2 wc find 



+i 



du 



VT^u 2 



K (VX) = tt I (a) K (b) ; 6>a>0 , 



(A9) 



+1 



du 



vT 



IoWX)=nI (a)I (b) 



(A10) 



du 



K {VX) = K (a)K (b) 



a, b > 



(All) 



du 



+i 



V^ 2 



I (VX)=I (a)K (b)+Ko(a)I (b) ; a,6>0 . 



(A12) 



Obviously, the first two integrals correspond to an integration over real angles ip, while the latter two integrals require 
for ip- values lying on the imaginary axis. The first three integrals are tabulated in mathematical tables like J27| ], but 
the integral ( A.12) seems to be a new result. Thus we shall shortly comment on a proof of this formula, which, by the 
way, also holds for the other integrals (A9) to (All): The idea is to show that the integral on the left hand side solves 
the differential equation of the modified Bcsscl functions with index zero with respect to a and b. After some partial 
integrations it turns out that this is ind eed the case; however, we thereby gather a boundary term, which exclusively 
vanishes for the boundaries chosen in (A9) to ( A12 ). For theses choices we then know that the integrals must be 
equal to a linear combination of the four possible products of Iq and Kq at argument a and b. The coefficients of 
these possible four contributions can then be determined by evaluating the integrals at a = b, and, if necessary, by an 
additional investigatio n o f the limi ts a — > and b — > 0. 

With the integrals (|A|) to ( |A12| ) we may now perform the </?-integral in ( |A8| ) and get a result of the form 



drj e^'Z^{a)Z^{b) 



(2), 



(A13) 



with Zq 1 ' and Z^' being modified Bessel functions with index zero. However, not any combination of Z^' and Z { 
with 77-borders at ±1 or ±00 will lead to a solution of the three equations 



r{2) 



,(2) 



d q 8 r - dp 



T q d r - O r d q 



a 2 - a 2 



(A14) 
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which necessarily should be the case for the limit A — > of the Chern- Sim ons integral because of ( 3.18 ). The reason 
for this is that, due to the complexification of the spherical coordinates ( |6.5| ), we might have integrated effectively over 
a three dimensional manifold S 3 , which produces boundary terms for partial integrations. For exam ple, there is no 
reason, why a surface described by re]l,+oo[, ipeiTl, rje] — 1,+1[ should lead to an integral solving ( |A14| ), although 
this integration contour leads to a finite integral. It will be a nice explicit check of solvability to investigate, for which 



choic es of Z^ and Z^ 1 and for which ^-borders ^>' cs according to (A.13) becomes a solution to the three equations 



7 (2) 



(Mi) 



While the operator T-$ annihilates the integrand of (A13) for arbitrary values of 77, and therefore ^' cs for arbitrary 
borders of the ^-integral, we have to perform suitable partial integrations to show the vani shing of T\ ^' cs and T2 ^' C s- 
Consequently, we get some boundary terms, which will not vanish for all possible forms of ( A13 ) . A d etailed discussion 
of these boundary terms reveals that there are exactly four possiblities to arrive at solutions of ( A14). For the adopted 
ordering of the three a v these solutions are given by the following integrals: 



d V e^I (a)K (\b\) 



1 it / &r,e°^I {a)I {b) 



(A15) 



+1 



123 



d V e^K (\a\)Io(b) , % 123 = in / dr, I (a) I (b) 



(A16) 



As an example, we want to discuss the action of the operator T\ on the wavefunction V E ,/ 12 , which can be written in 
the form 



J 1 ^ CS o 



' 1 " ' ' r ' ' 1 C 3 [h(a) I (b) + 7 (o) h(b)} - °2 lh(a) h(b) + J (o) I a (b)f 



(A17) 



That this is indeed equal to zero is easily seen for r\ = — 1, and, after some asymptotic expansions of the Bessel 
functions, also for r, — > —00. Similar results hold for the other solutions, where one of the two Bessel functions might 
be substituted by Kq or —K\. In such cases the vanishing of the boundary terms at r) = ±1 becomes a nontrivial 
question, because the K v - Bessel function might become singular at these points. A detailed investigation shows that 
there are only the four possibiblities given in (A15) and (A16) to avoid contributions from these boundaries. 

The question whic h rem ains t o be answered at this point is which integration contours for the (^-integral give rise 
to the integr ands in ( A15 ) and (A16). For the wavefunctions and 5 , o 12 3 the answer is easily found with help of 
the integral (A10), and we find the explicit representations: 



2tt 



*'l2 



lim 

2 7T A^O 



dr] / dtp 



: dr 



exp 



K{1,R) 



(A18) 



+1 



(112.", 



= lim 

2tt a^o 



dr] / dip 



K(1,R) 



exp 



(A19) 



where the Radius R of the r-integral has to be chose n with R < 2. However, for the other two wavefunctions the 
answer is more difficult, because, according to ( |A9| ) to ( A12 ), there is no easy possibility to generate their integrands 
uniformly for all 77-values of their integration regime. So we have to cut the jy-integral artificially at those 77-values, 
for which we get a = ±6, which actually happens at r\ = — crg/ci < —1 and r) = —<J\j<J2 > — 1. Suitable linear 
combinations of the results (A9) and (A12) then reveal the following integral representations: 
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dry J dip 
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123 2tt a^o 



+ 1 2 71 +00 — 2 7T +00 —<^l/o'2 +1 OO 

drj J dtp J —2 y drj J dtp J — J dr/ J dtp 
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X(1,_R) 



rdr 



exp 



(A21) 



Considering the representations (A18) to (A21) it is now possible to read off the integration manifolds E 3 which 
lead to solution s of t he Wheeler-DeWitt equation even for a non vanishing cosmological constant A. Although the 
representations (A20) and ( A21) look very strange because of the 77-junctions at —o%l<J\ and — 01/(72, they lead indeed 
to differentiable solutions. Since this is a non trivial claim, it will be proven in appendix C, where we show that the 
integrand of the outer 77-integral is a continous and differentiable function at the junction points, even for A ^ 0. 



APPENDIX B: ALTERNATIVE INTEGRAL REPRESENTATIONS FOR THE VACUUM SOLUTIONS 



We will now bring the one dimensional integral representations found in appendix A into a new, unified form, which 
will display nicely the symmetries of these solution s. Th e motivation of this new representation arises from the fact 
that one of the 77-integrals, namely ^0123 given in (A15), turns out to be a convolution integral, and may therefore 
be simplified with the aid of Laplace's convolution theorem. To become more precise, let us write 



1 

* / 0123 = 27T< J dr e - CT3T / [( ( T 1 -a 2 )T]e CT3 ( 1 - T )/o[(CTi+^)(l-r)] , 



(Bl) 



where we have substituted t = ^p- in (A16). If we denote 



g(a,f3;t) :=e at Io[Pt} , 



(B2) 



we may also write 



*0123 = 2lri 9(-<73,<T2 - (TV-, t) -kg (<T 3 ,(T2 + VV,t) | t=J 



(B3) 



i.e. (Bl) is a convolution integral, evaluated at t = 1. Let us recall Laplace's convolution theorem, which states 



C[h*h]{s) = C[h]{s)-C[h]{s) , with £[/](*):= J dtf(t) 





(B4) 



Using the inverse Laplace transform 



S + 2 OO 



f(t) = ^l J dse st jC[f](s) , 



(B5) 



where s has to lie right to all singularities of the integrand in (B5), Laplace's convolution theorem can also be written 
in the form 



f 



S+2 OO 
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With aid of the Laplace transform of g defined in (B2) 



C[g](s) = J dte^t I [pt] = [{u- s) 2 - [3 2 r h , 




(B7) 
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we arrive at once at the identity 



s+i oc 



Hi' - 

*0123 — 



J , — n 



(B8) 



if we employ (B6) with t = 1. Here we have introduced the four quantities s v via ( 5.16| ), which are ordered according 
to si < 32 < S3 < so for ci < 02 < 03. Since the quantity s occuring in ( JB8| ) must be placed right to all the s u , it 
must be chosen as s > so- With ( |B§| ) we have reached a very nice representation for x f , i23j which immediately shows 
up the symmetry with respect to arbitrary permutations of the s„: As the integrand obviously obeys this symmetry, 
and the integration curve, which is placed right to all the s v can be chosen to be the same after two of the s„ have 
been permuted, the integral must be invariant under such permutations. Translated to the <r p , we firstly have a 
symmetry under permutations of the a p themselves. Secondly, the additional symmetry under permutations of the 
form so <-> s p implies that the wavefunction ^0123 i s invariant under a reflexion <j\ — > —a\, <y% — > — (72,(73 — > (73 and 
cyclic permutations thereof. Thus we directly recover the symmetry properties of \E / q 123 , which we already claimed in 
section VI. 



Regarding the res ult (B8) in more detail, one may now ask, if there could be further integration curves, apart from 
the one used in (BE), leading also to vacuum solutions. Investigating the action of the operators T v defined in ( A14 ) 
we find the following result 
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and similar results hold for the action of the operators T\ and Ti on the s-integral. Thus we conclude that all four 
integration curves shown in figure [l] lead to solutions of T v ^' — 0, because the boundary term (|B9| ) vanishes at the 
end points of theses curves. However, up to now it remains unclear, how these vacuum solutions are connected to 
those which were found in appendix A. In the following, we shall be interested in a direct transformation between the 
two representations, which is rather tricky and therefore should be presented in detail at least for one of the remaining 
solutions. In our calculations, we will need the following two integrals: 



s-\-i 00 



V(a - s) 2 - /3 2 



s > a ± (3 



(BIO) 



J dte at K {P\t\] = -=f= 



^p 2 - a 2 

Let us now choose as an example the transformation of the wavefunction 
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To allow for the interchange of the r'- and the s-integral performed in the last two lines, we further have to make 
sure that I > s\, otherwise the r'-intcgral in the last line would not exist. If we further require s < S2, the r'-integral 
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will exist even for an extension of the r'-integral to — oo. Let us consider the additional contribution which we would 
get in case of this extension: 



s+i oo +1 

J dse s [(s-s )(s-s 3 )H J dr'e-^ +s >' K [(a 2 -ai)\r'\ 



+1 



s+i oo 



= / dT , e-—'K Q [{a 2 -a l )\T , \] / d S e^ 1 " r '> [{s - s )(s - s 3 )] ' = . 



(B13) 



Here the last s-integral vanishes, because, without meeting any singularity of the integrand, the integration contour 
can be deformed to the negative real axis, where the s-integran d van ishes exponentially. Thus, we still have a 
represent ation of ^>[, if we extend the r'-integral in the last line of ( B12 ) to the whole real axis. After employing the 
formula ( B 1 1 ) we arrive at 



1 



s-\-i oo 



ds [(s - s Q )(s - S3)] 2 [(s 2 - s)(s - si)] 



s+z 00 



/ dse s TJ ( s - s„) ^ , si < s < 



S-2 



(B14) 



Similar calculations reveal the remaining two identities 



s-H 00 



S2 < S < S3 , 



(B15) 



s+z oc 



*ias= / dseS f[(s-^) * > 



S3 < s < s 



(B16) 



It is now easy to comment on the symmetries of the wavefunctions ( B14 ) to (B16): While ^' 123 obviously is symmetric 
under permutations of the s p , and therefore under permutations of the a p , it is not symmetric under a permutation 
s p <-» so as the wavefunction ^0123- The vacuum state ^' 12 ^ s 0Iu y symmetric under the permutation Si <-* s 2 , and, 
consequently, under a\ <-» 02 ; the solution ^ has no symmetry of this form at all. 

Finally, we shall be interested in the limit Ti — > of our vacuum solutions, where we should remind the reader 
that the ^-dependence is hidden in the s„ via the a p , which depend on h according to (3.14). As an example, let 
us compute the limit Ti — ► of the wavefunction Vf^: First of all, the integration contour C\ can be deformed to the 
negative real axis, giving 



1 3 
*i = -2 f dse s f[(8-8„)-i 



(B17) 



Substituting s„ — Tis v and £ = Ti(s\ — s) we get 



*i = -2he s 



d£ 

7!' 



[(f~- Si + s 2 )(f - si + a )(£ - si + S3)]" 



(B18) 



In the limit Ti — > just an infinitesimal neighbourhood around £ = contributes to the ^-integral, thus we may expand 



2G 



*i -2fte Sl /" -^e-«> [(5i - s 2 )(s! - s )(*3 - Si)] * , (B19) 
o 

where we have substituted £ = in those parts of the integrand, which remain regular for h — > 0, £ — > 0, since we may 
choose e arbitrarily small. In the new variable £ = £/7i, we finally arrive at 

oo 

lim *i = , / ^ = -f- , . (B20) 

h -> Q V 2 (°"2 - ci)((T2 + 0-3X03 - ax) J V 2 ^/(o^ - <Ti)(a 2 + o 3 )(o 3 - a\) 

The limit % — > of the other vacuum solutions may be calculated analogou sly. As a result, we obtain the same 
asymptotic behavior as in the limit h — > 0, A — > 0, which we discussed in ( 5.10P of section V, i.e. the result is 
independent of the order in which these two limits are taken. 



APPENDIX C: CONTINUITY AND DIFFERENTIABILITY OF THE INTEGRAND ON THE 

ASYMMETRIC INTEGRATION SURFACES 



We finally want to show tha t the two integrals (A2C) and ( |A21 ) given in appendix A give analytical solutions of 
the Wheeler DeWitt equation ( 3.17 ) not only for A — > 0, but also for A ^ 0. To show this we have to prove that the 
effective integrand of the ^-integral is continous and diffcrentiable at the junctions rj = — 01/02 and r\ — — 02/01; the 
Wheeler DeWitt operator would otherwise produce boundary terms when acting on the wavefunctions and ^[23^ 
an d they could not be solutions. In the following, we will restrict ouselves to the solution \E ,/ 1 in the case A > defined 
in (A21), which can be written in the form: 



271-*; 



1 2-rr 

drj J dip 
b 



1 



exp 



-o^ja\ I 2tt 00 +i 00 

dr; I 2 J dip J + J dip 

-00 \ 1 — i oc 



K(l,R) 



r dr 
1 — r 2 



exp 



(CI) 



With the abbreviations P and Q defined via 

2-K ! 



P 



dp 



■ dr 



1 



exp 



+2 00 

dip 
-% 00 



K{1,R) 



exp 



(C2) 



we can rewrite (CI) as a piecewise defined ^-integral: 



-1 

27T*; = J dr] I 



P + 2Q , 7] < -02/0-1 
-P , rj > -<7 2 /<7i 



(C3) 



To show the continuity and differentiability of the integrand with respect to 77 at the 77-junction we then have to show 
the following two properties of P and Q: 



(P + Q)\ , =0 , — (P + Q)\ , =0 



(C4) 



To proceed in this direction let us try to evaluate the ^-integrals in P and Q. The ^-integral for P can be interchanged 
with the r-integral and is then peformed easily, leading to 



P = 2tt 



r dr 
1 — T' 



■Jo 



A , r 2 
— ao 

2 1 -r 2 



exp 



A 



(C5) 



with G' := 1 - r 2 + 2 K 3 rrj + -(«? + k 2 ) (1 - rj 2 ) + — 



4 1 -r 2 



2 KXK2V 

1 + r 



(C6) 
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The Q-integral is more complicated, because it is not allowed to interchange the r-and the ^-integral which have 
to be performed there: after such an exchange the (/j-integral would not exist. However, we are free to open the 
r- integration contour K(1,R) at r — > +00 , and the resulting integration curve can be deformed into a line integral 
ro — i R, where we can choose —1 < tq < 0. For this r-integration contour we may then interchange the r- and the 
(^-integral, because the (^-integral exists for any r-value along the contour. Thus we can write 



+z 00 



20 = 



1 -r 2 



dip exp 



2 1 



r dr 
1 — r- 



K 



[A , r 2 ' 




~G'~ 


— ab - 

2 l-r 2 _ 


exp 


A 



(C7) 



ro+i 00 



ro+i 00 



and in the last line the r-integral can be deformed back to the positive real axis, giving an integral around the cut for 
r > 1, which is generated by the -ftTo-Bessel function. Let us define this cut integral explicitcly in terms of the contour 
Cent via 



C C ut ■= K(l, R)® {r - iO 1 1 + R<r < 00} {r + i | 1 + R<r < 00} 



(C8) 



Obviously, C cu t is a superposition of a pole integral and two line integrals along the cut, where one of these line 
integrals is performed in the upper half plane, while the other one has to be evaluated in the lower half plane. The 
integral over the integration curve (CS) is independent of the Radius R, as long as this ranges between R = 2 and 
R = 0. Using the formulas 



K Q (-t±iO) = TiT-I (t) + K (t) , <>0 



(C9) 



we can evaluate the two line integrals which are contained in (C8), where it turns out that the l^Q-contributions 
cancel, while we get the io-contributions twice. Thus we can rewrite (C7) as 



2Q = -4tt 



1 -r 2 



l+R 



A , r 2 
— ab 

2 1 - r 2 



exp 



A 



2i 



1 dr 



K(l,R) 



1-r 2 



[A , r 2 - 




~G'~ 


[2 ah l-r\ 


exp 


A 



(CIO) 



If we now consider the limit R — > in ( C1C| ), the first integral takes the form of the P-integral, cf. (C5), and we get 

r 2 dr 



2Q 



-2P + 2i lim 



Jk(i,r) I-?" 



■K 



A , r 2 
— ab ^ 

2 1-r 2 



exp 



A 



(Cll) 



Comparing this result with the desired claim (C4), all that remains to be shown is that the pole integral in (Cll) and 
its ry-derivative vanish in the limit R — > 0. In this limit R — ► the only contribution to the integral arises from an 
infinitesimal region around r = 1, so we can expand the iio-Bessel function for large arguments; furthermore, parts 
of the integrand, which remain regular for r — > 1, may be substituted by their values taken at r = 1. We then find 
for the pole integral 



K(i,R) 



1 — 7" 



K 



A , r 2 
2 ab \— V 2 



exp 



A 



R^O 

(X 



dr 



K{\,R) y/l-r 



exp 



A (a-bf 
8 1-r 



(C12) 



where we have only taken into account the r-dependent parts of the integral, while a cr p -dependent prefactor has been 
omitted. If we now consider the right hand side of eq. ( |C12 ) at the r/-junction r\ = —cr^/cr-i, a and b become coincident. 
Thus the essential singularity of the integrand disappears, and all that remains is an integrable square-root singularity, 
for which the p ole in tegral vanishes in the limit R — > 0. Similar arguments sho w th at also the 77-derivative of the 
pole integral in (Cll) vanishes for R — > 0, so, after all, we have proven our claim ( |C4| ) and, therefore, the continuity 
and differentiability of the 77-integrand in ( |C1| ). Analgous calculations may be performed for the second non trivial 
integration surface Sf 2 3, where we have to deal with two different 77-junctions. 



(C13) 
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